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Abstract 

Triply heavy baryons are very interesting hadrons to be explored for they 
provide particular information about strong interactions, hadron structures, 
and weak decays of heavy baryons. We calculate the hadronic production cross 
sections of the £l ccc and the £l C cb, which are dominated by the gg fusion subpro- 
cesses containing 4362 and 1454 Feynman diagrams, respectively. A method 
for generating and calculating tree level Feynman diagrams automatically is 
developed to deal with complicated processes containing so many diagrams. 
Our results show that 10 4 — 10 5 events of triply heavy baryons can be accumu- 
lated for 10 fb~ l integrated luminosity at LHC. Signatures of the triply heavy 
baryons are pointed out, with emphasis on the decay modes f2 ccc — > Vt sss + 3tt + 
and il cc b — > Q sss + 3-7r + +7r~. We conclude that it is quite promising to discover 
them at LHC. 
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1 Introduction 



One interesting baryon to be discovered is the triply heavy baryon which is the 
color-singlet bound state consisting of three heavy c or b quarks by the strong inter- 
action. Especially, exploring this baryon will be very useful for understanding the 
three body static potential. Since the mass m of the heavy quark is so high, the 
relative motion of the quarks inside the triply heavy baryon is nonrelativistic and 
the typical velocity v of the heavy quark in the rest frame of the baryon is small. 
The three heavy quarks are bound tightly and the typical size of this baryon being 
order of 1/ mv is smaller than that of the conventional baryons with light flavors. The 
mass of the baryon M is the sum of the three heavy quarks' masses and the binding 
energy which is of order mv 2 . Similar with the heavy quarkonium system, in the limit 
of v <C 1, there are three distinct energy scales involved in the triply heavy baryon 
system, i.e., the mass of heavy quark m, the three-momentum of the heavy quark mv, 
and the off-shell energy of the heavy quark mv 2 . Nonrelativistic QCD effective theory 
(NRQCD) pp can be used to describe the triply heavy baryon system in a simpler way 
by reducing the number of energy scales, where physics effects at energy scale m are 
explicitly integrated out. Moreover, in the limits of both large N c and heavy quark 
mass, the mean field approximation can be applied to this system. Consequently, the 
wave function satisfies an ordinary differential equation [2]. Exploring this baryon 
certainly will enrich our knowledge about hadron structures and QCD interactions. 
In addition, the ground state can only decay via the weak interaction. Since the dy- 
namics is relatively simpler, the triply heavy baryon provides very interesting samples 
to study the weak decays of hadrons. In particular, both the identical fermion effects 
and the larger recoil effects caused by the spectators can be studied well in theory. It 
may provide some insight for clarifying some long-standing puzzles in weak decays of 
heavy hadrons. 

The masses, lifetimes, and other properties of these baryons have been studied over 
the past three decades SI El El El El O HH] • However, productions of these baryons 
are so difficult that they have not yet been discovered so far. They can be produced 
at high energy colliders by a direct production mechanism which can be described 
as follows. Triple heavy quark pairs are produced first at energy scale m or higher, 
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followed by the formation of the triply heavy baryons at energy scale mv . In the 
leading order approximation, three heavy quarks forming a color singlet state move 
together with smaller relative velocities v. Therefore, the production cross sections 
can be factored into short- distance coefficients and long-distance matrix elements. 
The short-distance coefficients describe the hard process effects of the triple heavy 
quark pairs production which can be calculated by perturbative QCD and expanded 
as a power series of a s at energy scale m or higher. The long-distance matrix elements 
describe the formation of the triply heavy baryon from point-like three heavy quarks 
with small relative velocity v. At hadron colliders, the triple heavy quark pairs can 
be produced via gg fusion and quark- ant iquark annihilation subprocess while at e + e~ 
colliders, it can be produced via e + e~ annihilation process. The production cross 
section for the triply- charmed baryon in the e + e~ collision has been calculated in Ref. 
[TT] and the predicted production rate turns out to be very small. LHC, however, as 
a very high energy and high luminosity hadron collider, provides a very good chance 
to discover this sort of baryons due to very large heavy quark pairs production rates. 

The authors in Refs. [12] and [32] calculated the hadronic production cross sections 
of the triply heavy baryons at LHC with their fragmentation functions. They derived 
fragmentation functions in the so-called diquark model of the heavy baryon and in 
perturbative calculation model, respectively. However, their obtained results do not 
represent the fragmentation functions for the production of triply heavy baryons 
predicted by QCD perturbation theory in the leading order. 

In [12], the author calculated the fragmentation function in the diquark model of 
the triply heavy baryons by treating two of the three heavy quarks as a point-like 
diquark. This is inappropriate for the reason that two heavy quarks can be treated 
as a point-like diquark particle only when the momentum of the gluon emitted or 
absorbed from the diquark is much smaller than the inverse power of the typical size 
of the diquark which is order of mv. In the production process, this condition cannot 
be satisfied since the momenta of those gluons are order of m or higher. Moreover, in 
his diquark description, the produced two free heavy antiquarks are forced to move 
in the same direction, which deviates from QCD description. 

In Ref. [13J, the authors calculated the fragmentation function using perturba- 
tion theory in Feynman gauge and including contributions only from two Feynman 
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diagrams. As shown in Ref. [2], the fragmentation function can only be correctly 
calculated in axial gauge. Moreover, there are seven Feynman diagrams contributing 
to the leading order fragmentation function which are shown in Fig.l in this paper. 

Even with the correct QCD leading order fragmentation function, one cannot ex- 
pect that it gives accurate description for the triply heavy baryons at small Pt- In 
Ref. [15], the authors did a comparative studied for the hadronic B c production by 
full QCD and by the fragmentation approximation. Their results show that the frag- 
mentation approximation is valid only when the transverse momentum of the heavy 
hadron, Pj-, is much larger than its mass M. However, in hadron collisions, the pro- 
duction cross sections are dominated in the smaller Pt region, typically comparable 
with M. Thus one can not expect that the fragmentation approximation gives accu- 
rate results. Instead, calculations of the production cross sections by full perturbative 
QCD are necessary and expected to give reliable predictions. 

In this paper, we study the production rates of the triply heavy baryons at LHC. 
We first calculate the short-distance coefficients in the leading order of ctf corre- 
sponding to the tree level contributions. Even at the tree level, there are a number 
of Feynman diagrams that need to be calculated and it turns out to be hard work. 
For instance, there are as many as 4362 Feynman diagrams for gg — > Q ccc + ccc and 
as many as 1454 Feynman diagrams for gg — > Q cc b + ccb. It is impossible to carry 
out the calculations by conventional calculation methods for cross sections. 

For a simpler production process e + e~ — > Q ccc + c c c, to simplify the calculation, 
the authors in Ref. [ 1 1 J took an approximation ignoring the heavy quark mass terms 
in the numerator of the quark propagator. The approximation significantly simplifies 
the calculations. However, as pointed out in Ref. [16], the approximation will lead to 
quite large errors with the same order contributions. 

To overcome the difficulty meeting here, we propose a method to evaluate Feyn- 
man diagrams by generating the diagrams automatically and use a direct amplitude 
calculation of the Feynman diagrams. The amplitudes are classified as a set of gauge 
invariant subsets in terms of independent color structures. Consequently, the code 
for numerical calculations is written in a significant compact form and numerical 
calculations of the amplitudes can be carried out. To ensure the correctness of the 
calculations, we verify various gauge invariance of the amplitude of each subset to 
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examine the code. 

Taking the long-distance matrix elements given in [5J, the cross sections of the 
subprocesses can be calculated. Incorporating with the parton model, we calculate 
the production cross sections of the triply heavy baryons at LHC. Our results show 
that 10 4 — 10 5 events of triply heavy baryons Q ccc and Q cc b can be accumulated for 
10 fb~ x integrated luminosity at LHC. We also present P^-distributions and rapidity 
distributions with various parameters. In this paper, we consider only the ground 
states with all orbital angular momenta vanishing. We also point out the signatures 
to reconstruct the triply heavy baryons, with Q ccc — > Q sss + 3n + and Q cc b — > fi sss + 
3n + + 7T _ being emphasized. We conclude that it is quite promising to discover them 
at LHC. 

The rest of the paper is organized as follows. In Sec 2, we present the calculation 
of the subprocesses of gg (qq) — >• Qq 1 q 2 q 3 QiQ2Q3- In Sec. 3, numerical results of the 
total cross sections and differential distributions are presented. Sec. 4 contributes to 
the signatures of these baryons and conclusions. 

2 Subprocess of gg (qq) -> Qq 1 q 2 q 3 QiQ2Q3 

In this section, we calculate the cross sections of the production subprocesses of 
Qq 1 q 2 q 3 at LHC in the leading order. There are a number of Feynman diagrams con- 
tributing to the short distance coefficients. We carry out the calculations by taking a 
strategy to generating the diagrams automatically and using a direct amplitude calcu- 
lations of the Feynman diagrams. Consequently, the code for numerical calculations 
is written in a very compact form. We present the calculations in detail below. 

2.1 Production mechanism of the triply heavy baryon 

The triply heavy baryon is the color-singlet bound state of three heavy quarks 
by the strong interaction. However, since the mass m of the heavy quark is much 
larger than Aqcd, the relative motion of the quarks inside the triply heavy baryon 
is nonrelativistic and the typical velocity v of the heavy quark in the rest frame of 
the baryon is small. The three heavy quarks are bound tightly and the typical size 
of this baryon being of order 1/mv is smaller than that of the conventional baryons 
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with light flavors which is of order 1/A.qcd- The mass of the state is the sum of the 
three heavy quark masses and the binding energy which is of order mv 2 . In the limit 
of i> <C 1, there are three distinct energy scales involved in the triply heavy system, 
i.e., the mass of the heavy quark m, the three-momentum of the heavy quark mv, and 
the off-shell energy of the heavy quark mv 2 . NRQCD [TJ can be used to describe the 
system conveniently, where physics effects at energy scale m are explicitly integrated 
out. 

The triply heavy baryon, £Iq 1 q 2 q 3 , being in a color singlet state implies that its 
color wavefunction must be of the form -^£^ 1 ^ 3 Qi£ 1 Q2&Q3t3i where (i— 1,2,3) are 
the color indices of the heavy quark Qi. The heavy quarks may be either the c 
quark or the b quark. The triply heavy baryons can be classified as two classes in 
terms of containing two or three identical heavy quarks. For those orbital angular 
momentum ground states, the exchange antisymmetry of the identical fermion implies 
that the triply heavy baryons with three identical heavy quarks can only be the spin- 
symmetrical states and only spin | states are allowed, while the triply heavy baryons 
with double identical heavy quarks can be either the spin-symmetrical state with spin 
3/2 or the other one with spin 1/2. We take Vt ccc , Q* cb and Q cc b to denote the triply 
heavy baryons consisting of the ccc quarks with spin 3/2, the ccb quarks with spin 
3/2 and spin 1/2, respectively. 

The mass of the triply heavy baryon is so high that it is very difficult to produce at 
normal high energy machines. At LHC, they can be produced by a direct production 
mechanism in which triple heavy quark pairs are produced first at energy scale m or 
higher by gg fusion or qq annihilation subprocesses, followed by the formation of the 
triply heavy baryons by combining the color-singlet triple heavy quarks moving in 
the same directions with small relative velocities at energy scale mv. In the heavy 
quark limit, t)<l, the production cross sections can be factored into the product of 
short-distance coefficients and the long-distance matrix elements. The short- distance 
coefficients describe the hard process effects of the triple heavy quark pairs production. 
They can be calculated by perturbative QCD and be expanded as a power series of a s 
at energy scale m or higher. The long-distance matrix elements describe the formation 
of the triply heavy baryons from point-like three heavy quarks with small relative 
velocities. In this paper, we consider only the contributions arising from the leading 
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order short distance coefficient and the leading order long-distance matrix element. 
There is only single long-distance matrix element. The factorization then holds at 
amplitude level which we illustrate below. Given the leading Fock state description 
presented in Appendix 16. 1[ by the standard perturbation theory, the amplitude of the 
subprocess gg (qq) -» ^QaQaQi^Os reads: 



Mag (qq) ->■ ^q 1 q 2 q 3 QiQ2Q 



3 J 



d 3 Vi d 3 V 2 V2M 



(2ir) 3 (2ir) 3 y/8rrii m 2 m 3 
^1^2) x M(gg (qq) -> (Q^Qst^Q^Qs) , (1) 

where M (M =mi+m2+rn 3 ) is the mass of the baryon Qq 1 q 2 q 3 , vrii is the mass of 
Qi (i=l,2,3), 4>(y \,V 2) is the wave function of the baryon ^Iq 1 q 2 q 3 in momentum 
space; V^i, V^ 2 and -(V^i + ^2) are the relative momenta among the three quarks 
in this baryon; d = 2, 3 is the number of the identical heavy quarks in the baryon 
arising from the identical fermion exchange antisymmetry and (QiQ2Q3) ( f' SzS> denotes 
the color-singlet state consisting of three heavy quarks with spin quantum number s 
and its third component sz- A heavy quark with a large mass moves in the triply 
heavy baryon with a small velocity v (in the baryon's rest frame). M(gg (qq) fusion) 
depending on \F X and V 2 describes the hard amplitude for producing the triple heavy 
quark pairs. The typical momenta are of order m or higher. Thus in the leading 
order (expansion in v) approximation, the 2 dependence in the matrix element 

M(gg (qq) fusion) can be neglected. The integration over variables V^i and \? 2 can 
be carried out. It follows that from equation (JTJ): 



M39 (99) ->■ ^QiOaOsQl^Qs) = 

[ M(gg (qq) (QiQ2Q3)r z> QiQ2Q 3 ) , (2) 



V2M #(0,0),, „, ,,, z . 



y/8m 1 m2 m 3 y/d). 

where \l/(0,0) is value of the space wave function of the baryon Qq 1 q 2 q 3 when the 
three quarks are all at the origin which describes the long-distance effects that happen 
at energy scale mv. In this way, the long-distance and short-distance effects are sep- 
arated at the amplitude level in the leading order approximation. Our task now is to 
calculate the short- distance amplitude M(gg (qq) — » (QiQ2Q3)i ,Sz ' 'O1O2O3)) which 
describes triple heavy-quark pair production with the triple quarks in a color-singlet 
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Figure 1: Seven basic diagrams for generating the Feynman diagrams of the triply 
heavy baryon production subprocess with different heavy flavors. 

and point-like {QiQ2Q^\ ,Sz ' > state. In calculating the amplitude, the momentum of 
each quark in the {QiQ2Qz)\ ,Sz ^ state satisfies Pq 1 : Pq 2 : Pq 3 = mi : m 2 : m 3 . 

In this paper, we consider only the triply heavy baryons production of the Q C cc 
0* c6 , and Q cc b states. Since the production rates of other ones such as flbbb, an d 
Q c bb are small, they are not included in this paper. 



2.2 Feynman diagram generating 

Since there are so many Feynman diagrams responsible for the amplitude of the 
subprocess, to carry out the calculation, it is crucial to generate the Feynman dia- 
grams automatically. We now propose a method to realize it by starting from the 
triply heavy baryons production via gg fusion subprocess. Notice that removing the 
external gluon lines and the corresponding interaction vertex for gg fusion subpro- 
cess, all Feynman diagrams reduce to seven basic diagrams as shown in FigJT] assuming 
three heavy quarks with different heavy flavors. They are nothing but the diagrams 
connecting three heavy quark lines in various ways. It means that all Feynman dia- 
grams for gg fusion subprocess can be obtained by inserting two external gluon lines 
into all possible positions of the 7 basic diagrams shown in Fig{D 

For the production process of Q ccc , the exchange anti-symmetry of the identical 
fermions means that the number of the basic Feynman diagrams is 7 x 3! = 42, while 
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for the production process of Q C cbi it means that the number of the basic Feynman 
diagrams is 7 x 2! = 14. 

We now show how to generate all Feynman diagrams for the triple heavy quark 
pairs production subprocess via gg fusion. Notice that every diagram contains at 
most single quadruple- gluon vertex. For convenience, all diagrams are divided into 
two classes in terms of whether they contain the quadruple-gluon vertex or not. 

We first generate Feynman diagrams without quadruple-gluon vertex by inserting 
the external gluon lines in between the gluon or quark lines of the basic diagrams 
one by one. There are 3 internal and 6 external lines for each of the basic Feynman 
diagram. The first external gluon can be inserted in between one of these 9 lines. 
Then there are 4 internal and 7 external lines. And the second gluon g can be 
inserted in between one of these 11 lines. Therefore, for the production process of the 
triply heavy baryon with three different flavors, the number of the diagrams without 
quadruple-gluon vertex is 7x9x 11= 693. 

We then generate the Feynman diagrams with one quadruple-gluon vertex. They 
are generated either from the last diagram in Fig{TJ where one gluon attaches to the 
triple-gluon vertex while the other one inserts into all other possible lines, or from all 
7 basic diagrams with the double external gluon lines inserting into the same point of 
any gluon line. These two ways generate 19 and 15 diagrams with quadruple-gluon 
vertex, respectively. Thus totally one has 693 + 34 = 727 Feynman diagrams. 

For the Q cc b (Q*cb) an d ^ccc production via the gg fusion subprocesses, there are 
727 x 2! = 1454 and 727 x 3! = 4362 Feynman diagrams, respectively. In this way, we 
generate all Feynman diagrams responsible for the tree level amplitude of the triply 
heavy baryon production via the gg fusion subprocesses. 

In a similar way, we now can generate all Feynman diagrams for the triple heavy 
quark pairs production via the quark-antiquark annihilation subprocesses. In the 
production subprocess of qq annihilating to three quark pairs with different flavors, 
Feynman diagrams can be classified as three types in terms of the number of gluon 
lines attached to the light quark lines. For the first type diagrams with one gluon 
attached to the light quark line, there are 7 x 9 = 63 Feynman diagrams without 
quadruple-gluon vertex and 1 Feynman diagram with one quadruple-gluon vertex. 
For the second type diagrams with double gluon lines attached to the light quark 
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lines, the Feynman diagrams can be generated by breaking up one of the gluon lines 
in the 7 basic diagrams and attaching to the light quark line. For this type, totally 
there are 2! x 15 = 30 Feynman diagrams. For the third type diagrams with triple 
gluon lines attached to the light quark lines, all gluon lines between heavy quarks 
are removed, and each of three gluon lines are attached to each of three heavy quark 
lines. Thus for this type, there are 3! x 1 = 6 diagrams. Adding the three types of 
Feynman diagrams together, totally there are 64 + 30 + 6 = 100 Feynman diagrams 
for the process of qq annihilation to three heavy quark pairs with different flavors. 
Thus for the Q cc b (Q* ccb ) and Q ccc production subprocesses via the quark- ant iquark 
annihilation, there are 100 x 2! = 200 Feynman diagrams and 100 x 3! = 600 ones, 
respectively. 

2.3 Color factors of the amplitude 

The large number of Feynman diagrams makes the calculation very complicated. 
To carry out the calculations we first need to account for the color factors. The 
amplitude can be classified into a set of gauge invariant subsets in terms of variant 
independent color factors. Denote & = 1,2, 3 , Xi — 1)2,3 (i = 1,2,3), a, b = 
1,2, •■ ■ ,8, as the color indices of three heavy quarks, three anti-heavy quarks, and 
two gluons, respectively. For the gg fusion production subprocess, there are totally 
12 independent color factors Ck (k = 1, 2, • • • , 12) listed as follows: 

-iii.ti-, k -Si4.;i ■ 

@1 = )?iX2%Xi%X3 > = (A A a ) ?lX2 % Xl % X3 , 

^3 = (A a A )£ 1X1 <% X2 %X3 > ^4 = (A A a )dxi%x2%x3 > 

c 5 = ^ (A a A )^ lX3 %xi%x2 > = (A A a ) ?lX3 % Xl % X2 , 
= (A a )axi(^ )6x2%X3 ) ^8 = — ^|-(A a )£ lX2 (A )6xi%xs > 

C9 = (A a )dX2(A )6X3%Xl ) ClO - (A a )6xa(^ )6X2%Xi > 

£?l?2i;3 £ Sl?2?3 

C11 = 77= (-^ixs^ )6xi%x2 > C12 = (A a )^ Xl (A )g 2X 3<% X2 • (3) 

The color factors given in Eq. ([3]) can be used as the independent bases of the color 
wave function of the amplitude of the gg fusion subprocess. With these color factors, 
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the amplitude of each diagram can be expressed as the sum of the products of the 
numerical values and those color factors. It yields: 

12 

M = d\J2c k T k , (4) 
fc=i 

where Tj, is the amplitude of the basic Feynman diagrams related to the color-factor 
C k . 

Denote q\ = 1, 2, 3 and ^2 = 1, 2, 3 the color indices of q and q for qq annihilation 
subprocess. There are totally 4 independent color factors D k (£=1,2,3,4) for the triple 
heavy quark pair production via qq annihilation subprocess listed as follows: 

^1 = /g "sixi "6x2 "6x3 > 

-^3 = /g "ciX3"«|l"?2Xl"$3X2 ' 

£ ?l€2?3 

-°4 = ^i?2%Xi%X2%X3 • ( 5 ) 

The color factors given in Eq. ([5} can be used as the independent bases of the color 
wave function of the amplitude of the qq annihilation subprocess. With these color 
factors, the amplitude of each diagram can be expressed as the sum of the products 
of the numerical values and those color factors. It follows that: 

4 

M = d\J2D k r' k , (6) 

fc=i 

where T' k is the amplitude of the basic Feynman diagrams related to the color-factor 
D k . 

For those diagrams without the quadruple-gluon vertex, the amplitude of each 
diagram can be expressed as single product of a numerical value and a color factor, 
each of which can be decomposed as the sum of the color factors listed above. For 
those diagrams with the quadruple-gluon vertex, the amplitude of each diagram can 
be expressed as a sum of three products of numerical values and color factors, again, 
each of which can be decomposed as the sum of the color factors listed above. 
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2.4 Simplification by identical Fermionic symmetry 

As discussed in Sec. 12. 2\ the permutation of the identical heavy quarks increases 
the number of Feynman diagrams by a factor 3! and 2!, for the production of color- 
singlet triple c quark state and double c quarks plus single b quark one, respectively. 
Consequently, there are 4362 and 1454 Feynman diagrams that need to be calculated 
for the amplitudes of Q ccc and Q cc b (fi* c b), respectively, via gg fusion subprocesses. For 
the quark-antiquark annihilation subprocesses, the number of the Feynman diagrams 
are 200 and 600 for the production amplitudes of the fl ccc and the Q cc b (fi^), respec- 
tively. Notice that the momenta of the identical heavy quarks in triply heavy baryons 
can be taken to be equal in calculating the amplitudes. This feature together with 
the identical fermion exchanging symmetry can be used to simplify the calculation. 
In fact, the contributions to the amplitude from those increased diagrams by each 
permutation of the identical heavy quarks are equal to the original diagrams. Thus 
one needs to calculate only those 727 diagrams generated via 7 basic diagrams for the 
gg fusion subprocesses, and those 100 diagrams for the qq annihilation subprocesses. 
Thus the calculations of the amplitudes can be simplified significantly. We emphasize 
here that it is available only for the amplitude of the leading order matrix elements. 

2.5 Automatic amplitude calculation 

With the simplification in the last subsection, the number of the Feynman dia- 
grams that need to be calculated is reduced to 727 and 100 for the production of the 
triply heavy baryons via the gluon-gluon fusion subprocess and the quark-antiquark 
annihilation subprocess, respectively. These are still very large numbers for doing nu- 
merical calculations. It would be quite lengthy and very tedious to code for numerical 
calculations if they are evaluated and expressed by the conventional method. Instead 
we develop a method to calculate Feynman diagrams closely related to the method for 
generating Feynman diagrams described in Sec. 2.2, incorporating a direct, numerical 
amplitude calculation method given in Ref. [17] as described as follows. 

Given the quantum numbers and momenta of the initial and final states in the 
amplitude, all momenta of the internal lines are fixed in a specific tree level Feynman 
diagram. The propagators and vertices in the Feynman diagram are described by 



12 



certain matrix or tensor while the spinor wavefunctions of the initial state and final 
state are described by the line and raw vectors, respectively!* Then the calculation of 
the amplitude of the corresponding Feynman diagram is nothing but the product of 
the matrices to a numerical number against Fermion lines. The method turns out to 
be efficient in numerical calculation of tree level amplitude. However, for subprocesses 
with too many Feynman diagrams like the triply heavy baryons production one, it is 
still difficult to write down the amplitude of all diagrams one by one. 

To overcome this issue, we develop an automatic calculation method, closely fol- 
lowing the method for generating Feynman diagrams described in Sec. 12.21 Notice 
that the Feynman diagrams differ from each other by the positions and the types of 
propagators and vertices. Each type is assigned to a fix number. A Feynman diagram 
is then represented by a table consisting of a set of ordered numbers, each of which 
corresponds to the type of the propagator or vertex in the diagram. With the quan- 
tum numbers and the momenta of each propagator or vertex, the Feynman diagram 
can be calculated by the multiplication of the matrices or tensors controlled by the 
number in the table. In this way, the program code for calculating the amplitudes can 
be written in a compact form and contributions from all diagrams can be calculated 
conveniently. In numerical simulation of each event, we assign the momenta, helicities 
of the light quarks, the heavy antiquarks, and the heavy baryon, and the polarization 
vectors of gluons randomly, assuming that they meet a set of conservation laws. We 
then use Vegas to simulate the events and carry out the phase space integrals. 

To ensure the correctness of the program code, we examine various gauge in- 
variance of the amplitude by Ward-Takahashi identities of the gg fusion subprocess. 
Namely, substituting one of two polarization vectors of gluons into its momentum, the 
total amplitude vanishes. It follows from Eq. Q that each amplitude r& (fc=l,2,...,12) 
vanishes under the substitution. To have additional test, in a similar way, we can 
also examine the triply heavy baryons production processes via double photons and 
'jg fusion subprocess, in which the amplitudes are expressed as the sum of different 
gauge independent subsets. Similarly, we can also have various test for the program 
code for the subprocess of the qq annihilation. 



^dn Appendix 16.21 we list the 4-dimension vectors describing the spinors of the quark and anti- 
quark, and the polarization vectors of the gluons. 
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2.6 The cross sections of the subprocess 

We have calculated the amplitude of the subprocesses gg (qq) to triple heavy 
quark pairs above. When the amplitude is expressed as Eq. (141) , using the formula 
given in Appendix 16.31 the squared amplitude for the gg fusion subprocess reads: 



+ — (r\r* + r\r* + r 2 r: + r 2 r* - r 3 r; - r 4 r* + r 7 r* + r 9 r* 
y 

32 

+ r n n 2 ) + — (r\r 7 - - r 2 r* + r 2 r; + r 3 r* - r 3 r^ - r 4 r; 
+ r 4 r* 2 - r 5 r; + r 5 r 12 + r 6 r? - r 6 r^ + r 7 r 10 + r?r* n + r s r; 
+ r s r* n + r 9 r* 2 + r w r* n ) - y (r^: + i\r* + - i\r* 2 
+ r 2 r* + r 2 r; + r 2 r* x - r 2 r* 2 - r 3 r* - r 3 r; + r 3 r* - r 4 r; 

rp* i "P "P* P 1 P 1 * i P 1 P 1 * P 1 P 1 * i P 1 P 1 * P 1 P 1 * P 1 P 1 * 
4 1 9 + 1 4 1 10 — 1 5 1 7 + 1 5 1 8 — 1 6 1 7 + 1 6 1 8 — 1 7 1 9 — 1 7 1 n 

4 

rp* p "P* P 1 P* P 1 P 1 * \ I / p p* "p p* p p* I p p* 

8 1 10 — 1 8-L 12 ~ 1 9 1 11 ~~ 1 10-L 12) + nl 1 l 1 8 " 1 I 1 9 ~ 1 2 1 7 + 1 2 1 1( 



+ r 3 r; - r 3 q 2 - r 4 r* + r 4 r* n - r 5 r* 10 + r 5 r* u + r 6 q - r 6 r* 2 )] 



where, Re means taking the real part of the number, and b means the summation 
of the color indices of the initial gluons and terminal heavy antiquarks. 

Similarly, with the amplitude given in Eq. (EJ), the squared amplitude for the qq 
annihilation subprocess reads: 



\m\ 2 = (d\) 2 3^r J r*-2- J Re(r' 1 r;* + r' 1 r;* + r;r;*) 



The differential cross section of the subprocess is then given by: 




k+l 



9 





dv{gg (qq) ->■ ^q 1 q 2 q 3 QiQ2Qs) 




E \ A te9 (qq) ^q iQ ,qsQiQ2Q3)\ 2 , (9) 



with, 



cffl 4 




d 3 P d 3 qi d 3 q 2 d 3 q 3 



(ki + k 2 - P - qi - q% - q 3 ) , 
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where s is the subprocess cm. energy squared, and the denominator N in 4 is 64 
and 9 for the gg fusion subprocess and the qq annihilation one, respectively, and the 
numerator 8 in -j| means the summation of the spins of the heavy antiquarks in the 
final states, which are given randomly in the code. The spins of the initial gluons 
(light quark and antiquark) are given randomly in the code too. 



3 Numerical results and conclusions 

Incorporating the parton model with the cross sections of the subprocesses, the 
cross sections of the process pp —> Qq 1 q 2 q 3 QiQ2Q3 + X in pp collisions can be calcu- 
lated. It follows that: 

a = J dx 1 dx 2 fg 1 (x 1 ,n F )f g2 (x 2 ,fJ,F) J da(gg ^ Q Qi q 2 q 3 QiQ 2 Q 3 , Vf) 

+ ^2 dx 1 dx 2 f q (x 1 ,n F )fq(x 2 ,HF) da(qq Qq 1 q 2 q 3 Q 1 Q 2 Q3,Hf) , 
g J J 

where / g (a;,/xp) and /Xp) (fq(x, /jlf)) are the parton distribution functions of the 

gluon and quark (antiquark) in the proton, respectively, /xp is the factorization energy 

scale, Ylq means summation of all light flavors of q. It happens that the predicted 

cross sections in the leading order are just the same for the production of the baryons 

consisting of identical fermions and non-identical ones. It is no longer valid beyond 

the leading order. 

Given the mass parameters m c = 1.5 GeV and m;, = 4.9 GeV, the values of 
|\I/(0, 0)|'s can be evaluated by the formula given in Ref. [6]. They are 0.0781 GeV 3 , 
0.0864 GeV 3 , for the Q ccc and Q cc b (fi* cfe ) states, respectively. 

Our numerical results show that the contributions to the total cross sections from 
the gg fusion subprocesses dominate as expected. In fact, contributions from the 
gg fusion subprocesses are 2~3 orders of magnitude higher than that from the qq 
annihilation subprocesses. Thus contributions to the cross sections from the q — q 
annihilation subprocesses are negligible. 

To carry our numerical calculations we take Cteq61 [18] parton distribution func- 
tions. Some uncertainties arise from the ambiguities in choosing the factorization en- 
ergy scale For comparison, here we take two different values of /xp, i.e., /xp=/x^/2 
and hf=^r, where jjL 2 R =P^+M . 
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Table 1: Predicted cross sections (in unit nb) of the triply heavy baryon production 
at LHC with v^S" = 7 TeV. Three typical Pt cuts are adopted. Pseudo-rapidity cuts 
|g| < 2.5 for CMS and ATLAS, and 1.9 < r] < 4.9 for LHCb are taken. 







LHC (CMS, ATLAS) 


LHCb 


- 


Vcut 

Prcut 


\v 


< 2.5 


1.9 < 7] < 4.9 








/V2 




/ifl/2 




OGeV 
5GeV 
lOGev 


0.0604 
0.00599 
2.6 E-4 


0.132 
0.0140 
6.3E-4 


0.0329 
0.00163 
4.8E-5 


0.0724 
0.00391 
1.21E-4 


0* 


OGeV 
5GeV 
lOGeV 


0.00151 
6.49E-4 
9.62E-5 


0.00351 
0.00152 
2.26E-4 


7.24E-4 
1.89E-4 
1.95E-5 


0.00172 
4.54E-4 
4.67E-5 


^ccb 


OGeV 
5GeV 
lOGeV 


4.89E-4 
2.43E-4 
4.49E-5 


0.00114 
5.67E-4 
1.05E-4 


2.15E-4 
6.86E-5 
0.894E-5 


5.09E-4 
1.65E-4 
2.13E-5 



We calculate the production cross sections of the baryons f2 ccc , Vt* cch and f2 cc b at 
LHC with total energy v^S* = 7 TeV and v^S* = 14 TeV and list in Table [1] and 
Table [2] with various Pt and pseudo-rapidity cuts, respectively. From the tables, we 
see that the predicted numerical results differ by around factor 2 with two different 
factorization scales ^f—^r/^ and ^f—^r- These uncertainties are expected to be 
reduced by including the next leading calculations which will be extremely difficult 
work. 

From Table [TJ we see that for 10/5" 1 integrated luminosity running at 7 TeV, for 
the production of f2 ccc , around (0.6 — 1.4) x 10 5 events at CMS and ATLAS can be 
accumulated with kinematic cuts Pt > 5 GeV and \rj\ < 2.5, while this number is 
around (2-4) x 10 4 with kinematic cuts P T > 5 GeV and 1.9 < 77 < 4.9 at LHCb. For 
the production of Q cc b, both fl cc b and VL* ccb needs to add together since the higher mass 
state will decay into the lower mass state by emitting a photon. After doing this, 
we see that for 10/6 -1 integrated luminosity running at 7 TeV, for the production 
of Q cc b, around (1 — 2) x 10 4 events at CMS and ATLAS can be accumulated with 
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Table 2: Predicted cross sections (in unit nb) of the triply heavy baryon production 
at LHC with v^S* = 14 TeV. Three typical Pt cuts are adopted. Pseudo-rapidity cuts 
M < 2.5 for CMS and ATLAS, and 1.9 < r] < 4.9 for LHCb are taken. 







LHC (CMS, ATLAS) 


LHCb 




Vcut 

PTcut 


M 


< 2.5 


1.9 < r] < 4.9 








/V2 




V.r/2 




OGeV 


0.113 


0.216 


0.0684 


0.135 




5GeV 


0.0123 


0.0258 


0.00412 


0.00906 




lOGeV 


0.000625 


0.00136 


0.000145 


0.000349 


0* 


OGeV 


0.00320 


0.00677 


0.00175 


0.00378 




5GeV 


0.00143 


0.00307 


0.000521 


0.00114 




lOGeV 


2.34E-4 


5.03E-4 


0.625E-4 


1.38E-4 


^ccb 


OGeV 


0.00105 


0.00222 


0.000527 


0.00115 




5GeV 


0.000544 


0.00117 


0.000190 


0.000419 




lOGeV 


0.000109 


0.000236 


0.289E-4 


0.639E-4 
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kinematic cuts Pt > 5 GeV and \r)\ < 2.5, while this number is around (2 — 6) x 10 3 
with kinematic cuts Pt > 5 GeV and 1.9 < i] < 4.9 at LHCb. 

From Table El we see that for lOOO/fr^ 1 integrated luminosity running at 14 TeV, 
for the production of Q ccc , around (0.6 — 1.4) x 10 6 events at CMS and ATLAS can 
be accumulated with kinematic cuts Pt > 10 GeV and \rj\ < 2.5, while this number 
is around (4 — 9) x 10 6 with kinematic cuts Pt > 5 GeV and 1.9 < rj < 4.9 at LHCb. 
For the production of Q cc b, again by adding both fl cc j, and f2* cb events together, we see 
that for 1000/fe -1 integrated luminosity running at 14 TeV, for the production of fl C cb, 
around (3 — 7) x 10 5 events at CMS and ATLAS can be accumulated with kinematic 
cuts Pt > 10 GeV and \r)\ < 2.5 while this number is around (0.7 — 1.5) x 10 6 with 
kinematic cuts Pt > 5 GeV and 1.9 < r\ < 4.9 at LHCb. 

Here we calculate only the contributions from the ground states. Excited states 
also can decay into the ground states. Including those contributions, the real pro- 
duction rates should be 2 — 3 times higher than ones listed in Table [1] and Table El 
Moreover, the cross sections for the anti-triply heavy baryons are exactly the same. 
This will double the accumulated triply heavy baryon events. From the tables, we 
see that the cross sections for LHC running at 14 TeV is about 2 times larger than 
that for LHC running at 7 TeV. 

We show the Py-distributions and the rapidity distributions of the productions of 
the triply heavy baryons in Figs. EH3 at LHC with various LHC running energies 
and in various detectors. The differential cross section drops fast with increasing Pt 
as expected. Some of the parts in the Pt distribution are not smooth due to Vegas 
simulation error which can be reduced by increasing the event numbers. 

It is instructive to look at the fragmentation description for the production of the 
triply heavy baryons. One expects that at large Pt the production is dominated by the 
fragmentation mechanism. The production cross sections are then factored into the 
product of the cross sections of the heavy quark pair production and the fragmentation 
function. As dimensionless quantities, the QCD leading order fragmentation function 
and the fragmentation probability must be proportional to en* |^(0, 0)| 2 /m 6 . Since 
the process is complicated, one cannot expect to have a simpler analytic expression. 

Comparing our predicted results for the Q ccc production with those fragmentation 
model descriptions presented in [121 [T3] , we see that our numerical result is around 2-3 
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orders of magnitude higher than that given in [12J and around 2 orders of magnitude 
smaller than that given in [T3] for > 10 GeV. These large discrepancies cannot be 
reconciled by varying parameters such as the heavy quark mass, factorization energy 
scale, type of the parton distribution function, and the rapidity cuts used in the 
calculations. 

To look for the reasons leading to so huge discrepancies between our result and 
those in [T2l 13], we first look at the nonperturbative parameter used in the calcula- 
tions. In |13j . the authors introduced a nonperturbative dimension-one parameter fg 
as used in the meson production when they calculated the amplitude of the production 
of the triply heavy baryon. From Eq.(5) in Ref. [13], the nonperturbative parameter 
fs they introduced can be related to our \l/ ccc (0,0) by |\I/ CCC (0, 0)| 2 ~ f\ ■ M^m^. 
In their calculation, they took j\ = 0.0625 GeV 2 corresponding to a value of the 
|^ccc(0, 0)| 2 in our paper being 0.458 GeV 6 , which is around two orders of magnitude 
higher than the value we used. From quark model calculations, this value is certainly 
over-estimated. This may explain the two orders of magnitude discrepancy between 
the results in [T3] and ours. 

In [12], the author introduced two nonperturbative parameters V'cc(O) and if) ; ccc (0) 
describing the wave function at origin of the double quarks forming the diquark and 
the quark and the diquark forming the baryon, respectively. His fragmentation func- 
tion was then proportional to lipcdfy'ipccdfyl 2 = 0.00345 GeV 6 , which was comparable 
to the nonperturbative parameter we used in this paper |\l/ ccc (0, 0)| 2 = 0.0061 GeV 6 . 
Thus the discrepancy between the results in [12] and ours cannot be attributed to the 
nonperturbative parameters used in the calculations. Notice that the fragmentation 
probability given in Eq.(8) in [T2] is proportional to a number A c he introduced. In 
the following equation after Eq.(8) in [12], the author showed that A c was a small 
number 4.2 x 10 -3 which arose from a cancelation between two very large numbers 
(Surprisingly it is a 10-digit cancelation!). This may give explanation for the dis- 
crepancy between our results and those presented in |12j . Physically, when they 
calculated the fragmentation function using the diquark model , probably only tiny 
contributions were accounted for. 

We also calculate the production cross sections for the triply heavy baryons at 
Tevatron. Our calculations show that there are around 2 thousands Q rrr events and 
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around 3 hundreds VL ccb events accumulated with 5 GeV P^-cut and 0.6 rapidity cuts 
for 2 fb' 1 integrated luminosity at Tevatron. Considering the detection difficulty, we 
think these numbers are not enough to discover the triply heavy baryons at Tevatron. 

4 Signatures and summary 

The particular signatures of the triply heavy baryons in the detectors can be 
used to reconstruct the triply heavy baryons. The ground states of the triply heavy 
baryons can decay only through the weak interaction. Naively, the decay width for 
Q ccc is triple the width of the free c quark while that for Q ccb is double the width of 
the free c quark plus the decay width of the free b quark. However, one expects that 
it receives large corrections both from the destructive interference between identical 
fermions and from large recoil effects. So the total decay widths may be deviated 
largely from the naive estimation. 

For fl ccc , a very interesting non-leptonic decay mode is that: 

Qccc 

I 

Qccs 

I 
Qcss 

I 

Qccc ^ Qsss ~t~ 3 7T 

Thus with this cascade decay mode, the Q ccc finally decays to Q sss + 3n + , which 
are stable charged particles in the detectors. Combining the tracks of the charged 
particles and the invariant mass, the triply heavy baryons can easily be reconstructed. 
For Q cc b, the similar but complicated interesting non-leptonic decay mode is that: 
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^ccb 







y 


V 








foes 


+ 7T + 




ccc 


+ 7T 






\ 




1 




^bss 


+ 7T + 




+ 7T _ 




+ 7T" 


1 








1 




^CSS 


+ 7T _ 




+ 7T + 




+ 7T 


1 




I 




1 






+ 7T+ 


^sss 


+ 7T + 




+ 7T" 



f2 cc & — > f2 sss + 37T+ + 7T~ 

Thus with these cascade decay modes, the Vl ccb final decay products are Vl sss + 
3 7T + + 7T _ , which are stable charged particles in the detectors. Again combining the 
tracks of the charged particles and the invariant mass, the triply heavy baryons can 
easily be reconstructed. 

We can also analyze the semileptonic decays of the f2 ccc , Vl ccb . The cascade decays 
are c — > s or b — >■ c — >■ s, with emitting leptons and neutrinos. This cascade decay 
feature with the decay vertex can also be used to identify the triply heavy baryons. 
Moreover, one can also combine the semileptonic decay with the nonleptonic decay 
to reconstruct the triply heavy baryons. 

In conclusion, our results show that a number of triply heavy baryons events can 
be accumulated at LHC. They can be reconstructed with their unique signatures in 
detectors. Triply heavy baryons are very interesting hadrons to be explored for they 
provide particular information about strong interaction, hadron structure, and weak 
decay of heavy baryons. They are still absent in the particle data booklet after the 
heavy quarkonium has been discovered over three decades. Our results show that it 
is quite promising to discover those triply heavy baryons in LHC experiments both 
for large number events and for their unique signatures in detectors. One may be 
waiting for an exciting time to discover them at LHC. 
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Figure 2: The predicted P<r-distributions of the triply heavy baryon production 
CMS and ATLAS, with v^S=7 TeV. 
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Figure 3: The predicted P^-distributions of the triply heavy baryon production 
LHCb with v^=7 TeV. 



23 




1E-6-I , , , , , , , , , , , , , , , , , , , , 

2 4 6 8 10 12 14 16 18 20 

P T (GeV) 



1E-3 q 




1E-6-| 1 1 1 1 1 1 1 1 , 1 1 1 1 1 1 1 , 1 1 1 

2 4 6 8 10 12 14 16 18 20 

P T (GeV) 



Figure 4: The predicted P<r-distributions of the triply heavy baryon production in 
CMS and ATLAS, with \/S=U TeV. 
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Figure 5: The predicted P^-distributions of the triply heavy baryon production 
LHCb, with \/S=U TeV. 
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Figure 6: The predicted rapidity distributions of the triply heavy baryon production 
at LHC with y/S=7 TeV. 
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Figure 7: The predicted rapidity distributions of the triply heavy baryon production 
at LHC with \/S=U TeV. 
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6 Appendix 



6.1 Nonrelativistic states 

The S— wave state of the nonrelativistic triply heavy baryon in its rest frame is 
constructed as: 



where, £ i; rji, and (Ei, pi) (i = 1,2,3) are the color, the spin, the mass, and the 
four-momentum of the heavy quark Qf, V\ and V2 are the relative momenta among 
the heavy quarks; (s, sz\r)i, 772, ^3) is the C-G coefficient; s and sz are the spin of the 
baryon and its third component, respectively; ip(vi, V 2) is the wave function of the 
baryon in the momentum space. 

The quark states are normalized by: 




(11) 



(f',Z',V,p\f,Z,V,P) = 5 ff 8^5 v > v (2Tr) 3 2E f 8®tf - p) , 



(12) 



where / and / are the heavy flavors. 
■00^1 j V^) is normalized by: 




(13) 



By the Eqs. (fl2]) and (ITB1) . we know the baryon is normalized as: 



(fi, s, s' z , P'\Q, s, s z , P) = (2n) 3 2E8 s ' s '5 Sz ' s 'z8 {3) (P' - P) . 



(14) 



The wave function at the origin is related to ip(Vi, V 2) by: 




(15) 



29 



6.2 Spinor and polarization vector 

For on-shell quark or antiquark with 4- momentum p^ = (p° , p l , p 2 , p 3 ) and mass 
m satisfying p 2 = m 2 , the spinor of the quark reads: 

uHp) = (h(\p\+P 3 ,P 1 + ip 2 ), f2(\p\+P 3 ,P 1 + ip 2 )) , 

u'Hp) = (f2(-p 1 + ip 2 ,\p\+p 3 ), M-p 1 + ip\\p\+p i )) , (16) 
and the spinor of the antiquark reads: 

vHp) = (-f2{-p 1 + ip 2 ,\p\+p 3 ) , fi(-p 1 + ip 2 ,\p\+P 3 )) , 
v'Hp) = (A(|p1+pV + *P 2 ) , -f2(\p\+p 3 ,p 1 + ip 2 )) , (17) 

where 



= yp°-|p| = yp° + |pi (18) 

v^bKbl+p 3 ) ' v^bKbl+p 3 ) ' 

Setting the beam line as the Z-axis, two physical polarization vectors of the gluon 

are: 

e+ = (0,1,0,0), e~ = (0,0,1,0). 

6.3 The A matrices and the Feynman rules 

In our calculations, the A matrices of SU (N c ) satisfy: 



[X a ,X b }=if abc X c , tr[X a ,X b ] = 5 ab , 

1 

~N~r. 



^(A a )y(A a )fei — 5a5jk — —SijSki , 



where N c = 3 in QCD. 

Feynman rules for those vertex are given as follows. 
Quark-antiquark-gluon vertex: 
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Tri- Vector vertex in order (123) with all momenta incoming to vertex: 

g_ 

Quadruple- Vector vertex: 




_|_jai03eja204e^Ati/x2^3/^4 _ ^/iiM^A t 2M3^ _|_ jaia^e ja%a^e ^i\[i2 gH3H4 _ g/J.i/j3 gH2V^j 
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